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The two-loop mechanism of Bjorken and Weinberg is used to constrain flavor 
changing neutral Higgs bosons. We calculate the complete set of two-loop diagrams 
for the rare decay fj, — > e+7 induced by such neutral Higgs bosons, for arbitrary Higgs 
and top masses. The analytic result is used to set limits on Higgs masses for some 
recent models with specific ansatz about the flavor changing couplings. For example, 
in the Cheng-Sher scenario of multi-Higgs doublet models, all neutral Higgs bosons 
possess flavor changing (/j <-> fj) couplings proportional to ^Jm^fnf. We find that 
the present limit on // — > cy implies that, in such scheme, these neutral Higgs bosons 
should be heavier than 200 GeV. 
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I. INTRODUCTION 



In most theories beyond the Standard Model, neutral Higgs boson couplings are typically 
flavor changing unless special arrangements, such as imposing discrete symmetries, are made 
to eliminate them. The mass of such Higgs bosons can be estimated from its potential 
contribution to the — K$ mass difference, Am^. If one assumes that the flavor changing 
sdH vertex has the same Yukawa coupling as that of the heaviest quark of the same type, 
that is, the 6-quark in this case, then ArriK implies that the Higgs mass should be at least 
150 TeV |]T[. This value is much larger than the electroweak symmetry breaking scale and 
immediately poses two potential problems: (1) the existence of an unnatural hierarchy in 
scale; (2) the Higgs sector would be strongly coupled and predictive power is lost. 

To avoid such pitfalls, two options are often exercised. The first is to impose some 
discrete symmetry to achieve what is called "natural flavor conservation (NFC)" |§, that 
is, to avoid tree level flavor changing neutral currents or couplings (FCNC). This is easily 
done by requiring that only one Higgs boson vacuum expectation value (VEV) contributes 
to each type of fermion mass ||. The second is to use some scheme such that tree level 
flavor changing neutral Higgs couplings Yij are naturally suppressed in low energy processes. 
The latter has received revived attention in the literature recently @||, in part because it 
may have interesting consequences at high energies, such as the decay |J of the top quark 
into the charm quark and a light neutral Higgs boson. 

For example, motivated by the Fritzsch ansatz of mass-mixing relations, Cheng and 
Sher pointed out [§J that low energy FCNC constraints may in fact be evaded in multi- 
Higgs doublet models without invoking the NFC condition. Let the contribution of the fc-th 

(k) 

Higgs doublet to the fermion mass matrix be , it is not unreasonable to assume that 
= x!>P ^Jm^fnj for every k, where is of order unity. Upon diagonalization, the 
fermion mass and mixing patterns can be roughly accounted for, but in general, neutral 
Higgs bosons in mass basis would all have flavor changing couplings. Their scheme can be 
summarized by assuming that the flavor changing couplings fi(F^ L P L + F^ R P R )fj H a of the 



2 



a'th neutral Higgs scalar H a have the following natural pattern 
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with A?- ' of order one. We shall refer to the lightest such scalar boson as H\. Since 
low energy constraints typically involve lower generation fermions, they are evaded by the 
associated tiny Yukawa couplings in Eq. ([[]). For example, taking A^ = ±A^ ~ 1, the 
ArriK constraint is weakened. If the lightest Higgs has pseudoscalar couplings, it is only 
required to be heavier than about 1 TeV, which is roughly the scale where the Higgs sector 
is expected to become strongly coupled. For scalar Higgs couplings, the bound is further 
lowered to roughly the symmetry breaking scale f8|. 

The Cheng-Sher scenario was not widely appreciated, and subsequent work || concen- 
trated on rare decays of r, B and D, as well as the utility of the [i — > c-f process as a 
constraint. It was recently realized, however, that the progressive nature of FCNC Higgs 
couplings could have interesting implications for very heavy quarks 0. In particular, the 
top quark may posses non-trivial couplings to a neutral Higgs boson H\ and a charm quark, 
of order m c t/v ~ ^m c m t /v. This may result in an appreciable branching fraction for a new 
channel t — > cH\ in the decay of the top quark. Alternatively, neutral Higgs bosons may 
have appreciable rates into tc type of final states. Because of the importance of top and 
Higgs physics, it is natural to ask whether we have exhausted low energy constraints on the 
Cheng-Sher scenario. 

In this paper, we give a careful analysis of the effect of these flavor changing neutral Higgs 
bosons on the \x — > e^ process up to the two-loop level. We parametrize our calculation 
in a way that is as model independent as possible so that our result can be applicable to 
any models with flavor changing neutral Higgs bosons. Our one-loop result differs from 
previous calculation || . Our two-loop study not only improves on previous rough estimates 
PJT0|j, but also uncovers some interesting characteristics that were overlooked before. In 
particular, the two-loop contribution can be larger than the one-loop result, as argued 
some time ago by Bjorken and Weinberg, and the heavy-Higgs-boson effect is not decoupled. 
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The contributions of the two-loop graphs that contain a heavy gauge boson loop diverge 
logarithmically as the Higgs mass is taken to infinity. However, due to a unitarity condition 
associated with flavor changing neutral Higgs bosons, if all such Higgs bosons are degenerate 
in mass, their contributions mutually cancel in such diagrams. This cancellation mechanism 
introduces model-dependence, reducing one's capability to make strict predictions. We 
present our results such that anyone can extract constraints on his favorite model. In this 
paper, We do not consider the charged Higgs boson, which gives an independent contribution 
with its own parameters. 



II. BJORKEN WEINBERG MECHANISM 

One flavor changing mode of particular interest is the celebrated fj, — > mode. The 
existing limit, at 4.9 x 10~ n ||11| , is one of the most impressive. It cannot occur at tree level, 
and it involves lepton number violation. In the Cheng-Sher scenario, the leading one-loop 
contribution has the neutral scalar and the r in the loop, with the photon radiated from 
the internal r line. The usual /i-/x-e sequence is much more suppressed. The r contribution 
(Fig. 1) to the 1-loop branching fraction of fi — > is 



BRl-loop(^ - *X) = I"" EKrKrTft^Tft + 



4vrm M 
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(2) 



where Mn a stands for the the mass of the a'th flavor changing Higgs boson H a . Clearly 
the contribution from the lightest neutral scalar boson dominates in general. This result 
holds if the flavor changing couplings are purely scalar (A^ L = A"- ) or purely pseudoscalar 
(A^ L = — A," R ). Also, our result is quite different from previous estimates Assuming 
lightest scalar dominance, we have 

B Rl _ loop („ e-,) = 5 x 1CT"|ALA^(HI^X)<(1 - 0.311n^!X) 2 , (3) 

which implies that, if one takes Ag T A^ r = 1 as in Ref. |J for the lightest scalar, any mass 
above 91 GeV is still phenomenologically allowed. This limit is more stringent than that 
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obtained in Ref. || because of the extra large \n(m 2 ./m 2 Hl ) term. The strong Mh 1 dependence 
in Eq.([|) means that the bound on Mh 1 will improve rather slowly with improvements on 
the BR(/i — > ej) limit. 

Here we wish to draw attention to an observation || made by Bjorken and Weinberg 15 
years ago, that certain two-loop graphs may in fact dominate over the one-loop contribution. 
The mechanism is as follows. Dipole transitions demand a chirality flip between the initial 
state and the final state of the fermion. For the one-loop graph involving virtual scalars, 
three chirality flips are involved: twice in the Yukawa couplings, and once in the fermion 
propagator. This fact is indeed just an accident at the one-loop level, but can be avoided at 
higher orders. Clearly, at the one-loop level, the \1-\1-e sequence is extremely suppressed, 
while even for the /x-r-e sequence discussed here, one pays the price of a suppression factor 
y/rnjri^m 2 / v 3 = O(10 -9 ). Going to two-loop order, one pays the typical price of g 2 /167r 2 , 
but if one could avoid two of the extra chirality flips, one may still gain enormously against 
the one-loop graph. In a set of two-loop graphs found by Bjorken and Weinberg, the virtual 
scalar boson couples only once to the lepton line, inducing the needed chirality flip. Through 
some heavy-particle (e.g. W or top) loop, the boson is then converted into two photons, 
one of which is reabsorbed by the lepton line. Assuming Mjj <C M^, Bjorken and Weinberg 
estimated the branching fraction due to the leading two-loop graph from the W contribution 
(Fig. 2a) to be 
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Here the neutral Higgs boson H a couples to the W boson at a relative strength cos 4> a with 
respect to that of the Higgs boson in the Standard Model. It is reasonable to assume that 
cos0 a is of order one, if the neutral Higgs boson originated from a Higgs multiplet that 
contributes significantly to SUl(2) breaking. Since this estimate works only when M# 
Myy, the present experimental limit implies that the Cheng-Sher scenario is probably not 
viable for light Higgs bosons, unless A" is significantly smaller than one, or the amplitudes 
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from different Higgs bosons cancel each other by accident, which could happen when the 
relevant Higgs bosons are degenerate in mass as discussed before. 

In the context of studying the electric dipole moment (edm) of the electron within neutral 
Higgs models of CP violation, Barr and Zee JT2| made independent observations that are 
analogous to that of Bjorken and Weinberg. Without knowing the work of Cheng and Sher, 



recently Barr [|Kj estimated the two-loop contributions to fi — > e^y for the case of very 
heavy Higgs bosons. Assuming only the W loop in the effective if 77 coupling and assuming 
Mfj 3> M^y, Barr estimated the branching fraction from the two-loop effect to be 



BR Barr (/, - ery) * - - - (cos^A^) 2 — ( 



ln Mr j 



(5) 
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Note that there is Mjj A suppression as in the one-loop case. Taking this result seriously 
and assuming cos0iA;L ~ 1 as before, one find that the present experimental limit requires 
Mh > 730 GeV in the Cheng-Sher scenario. 

Given the two estimates quoted above, where the lightest Higgs boson with flavor chang- 
ing couplings is either very light or very heavy compared to M^y, one may naturally be 
curious about the situation for Higgs masses in between. In this article, we report on a de- 
tailed calculation |T3[ of the complete set of two-loop diagrams where Higgs and top masses 
are kept arbitrary. 



III. COMPLETE TWO LOOP RESULTS 



The diagrams needed for the calculation of the transition dipole moment in our case are 



analogous to those for the electric dipole moments of quarks WM and electrons [12,15,16 



The most detailed depiction of these graphs can be found in Ref. JT6j, where one of the 
external electrons should be replaced by a muon. 

The two-loop graphs of interest can be classified into three sets, A, B and C, that are 
separately gauge invariant. Set A contains a heavy fermion loop. The fermion in the loop is 
most likely the heaviest one which is the top quark. They can be further classified into two 
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gauge invariant subsets. The first one involves an internal photon line in Fig. 2a, while the 
second is obtained by replacing the internal photon line with a Z boson line. According to 
Furry's theorem, only the vector coupling of Z boson contributes to the fermion loop. On 
the other hand, for both electron electric and magnetic dipole moments, the corresponding 
operators of the moments, a^ u and 0^75, are odd under charge conjugation, C, just as the 
vector coupling of Z. Therefore it is not hard to see that only the vector coupling of Z, 
not the axial one, to the external leptons can contribute to these operators. This argument 
can be carried over diagrammatically to the case of the transitional moments in the process 
H — > erf. Since the vector coupling of the Z boson is known to be relatively suppressed, the 
second group of diagrams can be ignored in the first approximation. 

Set B contains a W boson (and associated unphysical scalar) loop, and can be further 
divided into two gauge invariant subsets. For the first group, B/, the W boson loop induces 
a i?77 vertex as in Fig. 2. For the second group, B/j, the internal photon line is replaced 
by a Z line. Just like set A, set B n is suppressed compared to set Bj because of the small 
vector coupling of Z to charged leptons. In addition, if one assumes CP invariance, then 
since only the scalar components of the Higgs bosons couple to the W boson at the tree level, 
one expects that only the CP-even Higgs bosons will be relevant in this case. In general, 



from experience with the analysis of electric dipole moments in Refs. |T^JIJ-|l6[ , one expects 
set B to dominate over set A also. 

Set C involves graphs that have a different topological structure. They can be further 
divided into two gauge invariant groups Cj and C/j. They correspond to graphs without 



or with a Z boson line as in Figs. 3, 4 of Ref. respectively. Again, the second group 



is small compared to the first due to the small Z coupling. Numerical results of Ref. f!6 
indicate that, for the case of edm, the contribution of set C is in general much smaller than 
sets A and B. This conclusion should also be applicable to the transition dipole moment. 

We shall consider sets A and B first. For flavor changing leptonic processes, the internal 
gauge boson line can be either the photon or Z boson. However, if one is interested in flavor 
changing processes involving light quarks, a similar graph with both gauge bosons replaced 



by gluons (i.e. Hgg rather than i/77 or H^yZ) can also be important. 

The calculational strategy is to first calculate the one-loop effective vertex with one 
gauge boson, one photon and a neutral Higgs boson in the external lines. This has been 
done many times before ||17[| , and a recent calculation can be found in Ref. [|T^]. pseudoscalar 
amplitudes Higgs-gluon-gluon vertices Higgs boson are assumed to be on shell in Ref. flIB| , 
it is easy to extract the result with off-shell Higgs boson as long as one can tell which factor 
of Higgs mass comes from the loop momentum and which one is due to the vertex. The 
result of Ref. is consistent with the recent calculation of electron edm [16H, where the 
Higgs boson was kept explicitly off-shell, but only i/77 contribution was given. 

Here we shall concentrate on leptonic FCNC In that case the graphs with internal Z 
boson are suppressed relative to the ones with internal photon line by a factor of (1 — 
4 sin 2 6 l i4/)/4sin 2 6 Wl which is about 0.087 for sin 2 W = 0.23. Therefore one could ignore 
these contributions even though they can be easily incorporated into the analysis. 

We shall parametrize the relevant couplings as 



C 



m t 



^ t(A° t P L + A a tt *P R )t H a 



V ^ 6 KKuPl + K»Pr)H H a + gM w cos <j) a W + W~H a 



(6) 



Here v = (V2G F )~? ~ 246 GeV. If one imposes CP conservation then Im (A" f ) 2 = and 
Im (A^JA^J*) = 0. Note also that in case of CP conservation, cos0 a is nonzero only for 
those CP-even scalar Higgs bosons. (For the Higgs boson in Standard Model, A u = 1) 

To simplify long expressions, we define the reduced amplitude A, which is dimensionless, 
for the transition fi —>■ ej(e, k) as follows: 
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and the branching fraction is 
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Note that CP conservation will require Im (AlA* r ) = 0. For set A with the top-quark 
loop, the i/77 or H^Z vertices already contain one power of external photon momentum. 
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Therefore, we can set the virtual photon momentum and the Higgs boson momentum to be 
equal and the two-loop result can be easily produced. The i/77 contribution gives 

where z ta = rn^/M^. The chirality factors are defined as A5 = — 1 and A,? = 1. The scalar 
Hti Yukawa coupling Re A° t is associated with the following function, 

f(z) = -z dx \ ^ln^ ^. 10 

2 jo x(l — x) — z z 

The pseudoscalar coupling Im A° t is associated with 

#(z) = -z t dx- l - In x ^~ x \ (11) 

yK ' 2 Jo x(l-x)-z z K ' 



We have closely followed the notations of Ref . WW . 



If CP is invariant, the a'th Higgs boson, when it couples to the top quark, is either a 
scalar (Im A" t = 0), or a pseudoscalar (Re A" t = 0). However, we do not have relations 
between A^ L and A" except that they are relatively real. 

For the Z-mediated diagrams, 

a l,r HZy s _ (1-4 sin 2 6 w )(l - 4Q t sin 2 W ) 

^Vloop ^^ e + 7J 16sin 2 ^cos 2 ^ 

x3Q t J2 A^ L,fl 2[Re A a J(z ta) z tz ) - zA^Im A%g(z ta , z tz )] . (12) 

a 

Here f(x,y) = yf(x)/(y - x) + xf(y)/(x - y) and similarly g(x,y) = yg(x)/(y - x) + 
xg(y)/(x — y). We have also extended the previous definition to denote z t z = m^/Mf. Note 
that, in this Bjorken- Weinberg mechanism, there is only one power of light quark mass 
suppression ||, which has been explicitly written out in Eq.(|7|). 

To derive the contribution of the bosonic loops, we shall classify the graphs into two 
gauge invariant types. The first type does not depend on Higgs mass in their couplings 
while the second set does. As a result, the first set is power suppressed by the Higgs mass 
while the second set is logarithmically increasing when the Higgs mass becomes very large, 
which is a very intriguing situation. For Higgs mass larger than a certain value the second 
type dominates. We shall only present the combined contribution of the two types. 
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The H a WW vertex in Eq.fl6|) is parametrized as gMwg^ u cos0 a , where cos0 a is a Higgs 
mixing angle. For the Standard Model Higgs boson, cos</> = 1. Before we proceed with 
our results, it is important to point out a unitarity constraint on the flavor changing neu- 
tral couplings. One can always make linear combinations of the scalar doublet fields such 
that only one doublet is responsible for symmetry breaking. It is the scalar component 
of this doublet that couples to W boson pairs. Since this combination is also responsible 
for generating masses to the fermions, its Yukawa couplings should be automatically flavor 
conserving. Upon diagonalization of the Higgs boson mass matrix, all neutral Higgs bosons 
should in general possess flavor changing couplings. However, the above observation leads 
to a unitarity condition 

y>os0 a A* L ' fl = O (fori /J), (13) 

a 

which basically reflects the fact that the scalar doublets that mediate flavor violation must 
have zero vacuum expectation value at tree level. One important consequence is that, for 
the graphs in set B, terms that are independent of Higgs mass are cancelled away. 
For the Hjj case, one obtains the two-loop amplitude 

A w R \oo> -> e + 7 ) = - cos <f> a A^ R + 5g(z a ) + \g{z a ) + \h{z a )] , (14) 

with z a = M^/Mfj^ The function h(z) is defined as 

r/\ 2 9 f g{z)\ z r 1 dx r z . x(l — x)~\ 

h( z ) = z 2 —[^— L ) = - - -H - -ln-^ -. (15) 

oz\ z J 2 Jo z — x(l — x) L z — x(l — x) z 

It is straightforward to see that A L,R oc Ag M L,fi for this H^-loop amplitude and other 

purely bosonic loop contributions, unlike the situation for the t-loop with the pseudoscalar 

coupling ImA" t . Therefore we shall drop the chirality label for amplitudes from sets B and 

C for brevity. Note that if CP is conserved, cos <p a = unless the H a is a scalar. Therefore 

the flavor changing pseudoscalar Higgs boson does not have contribution of this class in such 

case. 

It is useful to know the shapes of these functions /, g and h. Numerically, /(l) is about 
0.8, while g(l) is about 1.2. The general z dependence of these functions are given in Fig. 3. 
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Unless z is very small or very large, these functions are of order unity. For very large or 
very small z p^JI^ JiT 



f(z > 1) ~ §lnz + #(z> 1) ~ §lnz + l, h(z > 1) |(W + 1), 

/(z < 1) ~ f(hiz) 2 , 1) ~ §(lnz) 2 , h(z <£l) ~ zlnz. (16) 

For the large z asymptotic forms, we obtain Eq.(£|) in the light-Higgs limit from Eq.(|l4l). 
It is tempting to use also Eq.(14) to find the heavy-Higgs 2; < 1 limit, which will produce 



the estimate Eq.(y) given in ref. ||10|| . However, this estimate clearly overlooks other non- 
decoupling contributions in Fig. 2c, d that we will discuss. 
For the HZ'-f case, one has 

x |(5- tan 2 9w)f(z a ,Zz) + |(7 - 3 tan 2 9 w )g(z a , z z ) 

+\g(z a ) + \h{z a )] . (17) 

with zz = M^y/M^. The suppression factor of (1 — 4 sin 2 9 w )/4sm 2 W comes from the 
vector part of the eeZ coupling. The W-loop contribution of HZ'y is about 10% of that of 
and they have the same sign. 

If we assume that CP is conserved and the Higgs boson is a CP-odd pseudoscalar, it does 
not couple to the W boson and there will be no contribution from set B and C. However, 
there will be ImA^ contributions in Eq.(||D from set A. 

The type of diagrams that involve the Higgs mass squared in the coupling are shown 
in Fig. 2b,c,d, which are Fig. 2k,l,m in Ref. 0] respectively. They contain the coupling of 
the physical Higgs boson H to the unphysical Higgs pair G + G~ . An important exception 
is the contribution related to Fig. 2b. This diagram has been grouped into Fig. 2a with the 
bosonic inner loop because they are gauge related. The contribution proportional the Higgs 
mass squared, from the HG + G~ vertex, can be combined with part of Fig. 2a to form a 
vertex that is proportional to the inverse Higgs propagator which then cancels with the Higgs 
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propagator in the outer loop. Thus the resulting contribution for each H a is independent of 
M Ha and therefore cancels each other completely because of the unitarity condition Eq. (|T3|). 
The contribution from Fig. 2c, d gives 



^G-IOOP^ - e + 7 ) = - E COS <PaK,— \f(Za) ~ 9 (z a 



and 



< Z loop(^ e + ^) 



2z a 



(18) 



1-4 sin 2 9 



w . 



8 sin 9 W 



tan 2 W ) 



xEcos0 a A: — 

a Lz a 



f{z a , Z Z ) ~ g{z a , Z z) 



(19) 



Note that /(l) — g(l) = —0.4, while for small z, which corresponds to the case of very 
large Higgs mass, f(z) — g(z) ~ z(lnz + 2). This leads to the peculiar situation where 
the contribution increases logarithmically with the Higgs boson mass. The coefficients are 
small enough that these contributions are not so significant as compared to the Vy-loop 
contribution discussed earlier, except for the case of very heavy Higgs boson. Of course, one 
may not trust the perturbative estimate if the Higgs mass becomes too heavy and the Higgs 
self-coupling becomes nonperturbative. 

The contribution of the two-loop graphs in set C can be easily translated from the 



calculation of Ref. JT6|. The result is 

1 



A c (n -> e + 7) 



4 sin 9 W - 



(z a ) + Df ) (z a ) + (z Za ) + Df h ) (z Za ) + D&> (z Za ) 



(20) 



where the functions D^ 3a ' b ' c ' d ' e \z) and D^ a,c \z) are given in Appendix B of Ref. [16], zz a = 
M z/ M H a and D( e b \ z ) = 4: sin 2 6 W tun 2 6 w D^ c \z). Note that the terms with functions 
D (3a,b,c,d,e)( ZH } belong to the group Ci, while the rest belong to the second group Cjj. As 
commented earlier the second group is suppressed relative to the first group. The reason 
one can easily translate the calculation of edm from these graphs into contributions to the 
transitional magnetic moment is because the Higgs line is always attached to one of the 
external fermion lines, and because the Higgs boson only has scalar couplings to gauge 
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bosons at tree level. Therefore, in the case of edm the Higgs coupling to fermions is always 
pseudoscalar, while its coupling for transitional magnetic moment is always scalar. 



Just as the calculation of edm in Ref. fig] , the contribution of the graphs in set C to 
the transitional magnetic moment is also small. We therefore do not include them in our 
numerical analysis. 



IV. DISCUSSION AND CONCLUSION 

In our numerical analysis, we shall ignore CP violation and take Af t to be real. We 
also assume that A^ = = A^, i.e., scalar Higgs couplings. Under this condition, the 
reduced amplitudes are simplified A L = A R = A. Comparison with experiment is given 
in Fig. 4. The result in general depends on three parameters, in addition to the unknown 
top and Higgs masses. They are the parameters A a and A eA1 (Eq.(H)) which parametrize 
the ttH and e\iH couplings, and cos</> a , which parametrizes the H a WW coupling. We set 
these parameters to one in our figures as a reference point. At this moment, we pretend 
that the contributions from different Higgs bosons do not strongly cancel each other. The 
results due to the contributions from a single Higgs boson are shown in Fig. 4. Numerically 
the contribution from the top-quark loop is generally smaller than that from the iy-boson 
loop. Also, the contributions due to the Z boson can be ignored although we have included 
them in our numerical analysis. 

It is instructive to look at the numerical results at the amplitude level. In Fig. 5, we show 
the reduced amplitudes A's from various sources. The error bar of the data point indicates 
the experimental bound on /i — > ej, \A\ < 1. The data point is used to guide the reader's 
eyes. It is purposely located at the lower bound of the Higgs mass at 91 GeV from the 
one-loop result. Because the two-loop W contribution does not vanish in the heavy-Higgs 
limit, one may need to sum up contributions from different Higgs bosons. For the case of 
the two-doublet model, the unitarity condition implies that the W amplitude is just the 
difference between those from the Higgs bosons at two separate masses. It is understood 
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that when the Higgs bosons are degenerate in mass, the WMoop contributions cancel each 
other. 

Some interesting features deserve special attention: 
(1) As Mu —> oo, the amplitude in Eq.(^) does not go to zero, unlike Barr's estimate in 
Eq.([5p. In fact it goes to infinity as lnM#. This non-decoupling behavior is curious but can 
be easily understood. In Feynman gauge, the non-decoupling graphs involve neutral Higgs 
coupling to unphysical charged scalars, which is proportional to Higgs mass squared. Such 
couplings are dictated by the gauge symmetry and its breaking. At tree level, the Higgs field 
that appears in the W + W~H coupling is the scalar component of the Higgs doublet that 
generated the symmetry breaking (the other three components are precisely the unphysical 
scalar bosons). This component is in general not a mass eigenstate of course. It can be 
expressed as the sum of the scalar components of the Higgs fields that participate in the 
breaking of SUl(2), with coefficients proportional to their respective VEV's. Therefore any 
Higgs boson that couples to the W pair has to originate from some weak multiplet that 
contributes significantly to the breaking of SUl(2). Hence the natural scale for such Higgs 
bosons should be the SUl(2) breaking scale, v. If their masses are artificially pushed much 
higher, say, by fine tuning, their physical consequences would not decouple. For similar 
reasons one can also conclude Jl9| that, in Standard Model, the Higgs contribution to g — 2 
of charged leptons would also not decouple. For very large Higgs mass, it should diverge as 
In Mjj. 



(2) In set B, the contribution from Eqs. fll8yi9D that gives rise to the non-decoupling behavior 



has a different sign compared to the other part from Eqs.( |T^|T7|) . At a low Higgs mass below 
200 GeV, the non-decoupling contribution are small. This is the region where Barr's estimate 
applies []nj. However, around 600 GeV a perfect cancellation occurs as shown in the solid 
curve Fig. 5. Even if the contribution from the top-quark loop is included, it will only shift 
the position of the amplitude zero by a small amount, depending on the value of m;. If the 
mass of the flavor changing neutral Higgs boson happens to lie in this region, BR(/i — > ej) 
would be very suppressed and further improvements of the experimental limits will not be 
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very constraining. Such cancellation behavior can potentially become a crucial issue in the 
future. 

(3) The present experimental limit requires Mh > 200 GeV under the same simplifying 
assumptions made earlier. This is a factor of at least 2 better than the one-loop limit. Our 
conservative limit is substantially lower than the bound 730 GeV from Eq.(^), where only 
some of the W-loop contributions are included. It turns out that other contributions from 
the non-decoupling term and the t-loop diagram reduce the bound substantially. 

(4) The result is only mildly sensitive to the top quark mass. However, there exist models 
in which the flavor changing neutral Higgs bosons have very small couplings to the W boson. 
That is, cos may be very small. In that case the top quark contribution dominates at two- 
loop and is still larger than the one-loop result. Conversely, there are some other models in 
which the Higgs boson that couples to leptons is different from the one that couples to up- 
type quarks. In that case, the parameter A tt would be zero and the top-loop contribution 
should be ignored. 

(5) If the Higgs boson couples also to the down type quarks as in the general Cheng-Sher 
scheme, the constraint from Kl — Ks mass difference in general would dominate over those 
from fi — > ej, although they are modulated by different Ay factors. 

(6) The limit on /x — > e^ of course will improve in the future. However, note that most 
of the severe low energy FCNC constraints such as K^-Ks mass difference, B°-B° mixing, 
Kl — > fie and \i — ► ej etc., originate from processes involving down-type quarks and charge 
leptons. FCNC constraints involving up-type quarks (e.g. D°-D° mixing) are rather weak. 
As pointed out in Ref. ||, it is in fact easy to avoid constraints from K, B and /i systems 
completely, by assuming NFC for down-type quarks and for charged leptons. Although 
the ansatz may seem a bit artificial, it does, however, permit tantalizing phenomenological 
consequences for the top quark ||, despite the depressed \i — > ej transitions. 

To conclude, we have derived the result for the two-loop contribution of flavor changing 
neutral Higgs bosons to the celebrated rare decay \x — > ej, for arbitrary Higgs and top masses. 
This is one of the rare situations in which higher order contributions actually dominate over 

15 



lower order ones. The numerical consequences depend, of course, on the model. For the 
generic case in the scheme of Cheng and Sher 0, the result shown in Fig. 4 improves the 
one-loop bound by more than a factor of two. The curious behavior of non-decoupling of 
very heavy Higgs boson effects at two-loop is emphasized. 
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FIGURES 

FIG. 1. A one-loop Feynman diagram for fi — > e + 7 through the r lepton as the intermediate 
fermion. 

FIG. 2. Feynman diagrams for ji — > e + 7. The generic inner-loop in (a) involves the t-quark, 
the IF boson and its ghost. For purely bosonic contributions, diagram (a) includes sea-gull graphs 
and other gauge related graphs except that we separate out those with vertices G + G~ H oc M\ 
in different diagrams (b), (c) and (d). We also do not show conjugate diagrams with lines of the 
neutral gauge boson and the Higgs boson exchanged. 

FIG. 4. NumMcifl^flniai^ a ^lM*kM fine 1 Me one-loop result 

via the intermediate r lepton. The two-loop result, assumming coming from one single neutral 
Higgs boson at M#, is given by the solid (dashed) curve for the case m t = 100 (200) GeV and 
cos 4, A" = 1, A tt = 1. 

FIG. Ti. Reduced amplitudes A's for the process \x — > e + 7. The dotted line is the one-loop 
result via the intermediate r lepton. The data point is located at the Mh lower bound due to the 
one-loop result. The error bar indicates the experimental bound on the reduced amplitude. The 
two-loop t contribution is given by the dashed (dash-dotted) curve for the case mt = 100 (200) 
GeV and A tt = 1. The two-loop W contribution, for the case cos </> a Ag^ — 1, is given by the solid 
curve, which does not vanish in the heavy-Higgs limit. 
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